Quantum dynamics of spatial decoherence of two atoms in a ring cavity 
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We study the spatial decoherence dynamics for the relative position of two atoms in a single- 
mode ring cavity. We find that the spatial decoherence of the two atoms depends strongly on 
their relative position. Taking into account the spatial degrees of freedom, we investigate the 
entanglement dynamics of the internal states of the two atoms. It is shown that the entanglement 
decays to almost zero in a finite time, and the disentanglement time depends on the width of the 
' wave packets describing the atomic spatial distribution. 
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Q . I. INTRODUCTION 
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' Superposition and its many-particle version — entanglement — are two basic featmes of quantum physics, distin- 
guishing the quantum world from the classical world. Because of its intriguing properties, quantum entanglement 
■ has attracted considerable attention as an important resource for quantum information processing [HQ- However, 
quantum coherence can be destroyed due to the physical system interacting with the environment, which has been 
recognized as a main obstacle to realizing quantum information processing. Hence a better understanding of the 
C ' mechanisms of quantum decoherence is not only crucial for the understanding of the quantum-classical transition 
(see, e.g., [3,-^^]), but also essential for the implementation of quantum information processing. 

In the last few years, theoretical studies in this context have involved a variety of systems (see, e.g., Refs. j6l-[l^). 
Moreover, experiments have also been done to demonstrate the dynamic process of decoherence as well as the collapse 
and revival of the quantum coherence (see, e.g., Refs. [l3l - [l5| ). In recent years, several physical systems have been 
' studied to learn more about environment-induced decoherence (see, e.g., Refs. 

. The influence of atomic spatial motion on quantum dynamics has been considered in different contexts (see, e.g., 
' Refs. [H, [l^). For a system of two cold atoms placed in a noisy vacuum field, the back-action of emitted photons 
, on the wave packet evolution about the relative position of the two cold atoms was discussed in Ref. [l^. It was 
' shown that the photon recoil resulting from the atomic spontaneous emission can induce the localization of the relative 
' , position of the two atoms, through the entanglement between the spatial motion of individual atoms and their emitted 
photons. 

, In contrast with previous works, here we consider the case where a pair of atoms simultaneously interact with a 
T-H . single- mode ring cavity. In fact, this kind of system is a typical system in which the cavity field acts as a data bus 
^ ■ inducing two qubits to be entangled (2TI . l22j . The system is efficiently used in many schemes for proposing how to 
• i-H realize quantum logic gates and teleportation in cavity QED [23| . In this paper, we will study the spatial decoherence 
of the atomic relative position and the disentanglement of the internal states of these two atoms induced by the 
^ , back-action of the photons emitted from these atoms. We find that the spatial decoherence of the two atoms depends 
■ - - ' strongly on their relative position. Our results show that the entanglement of the internal states decays to almost 
zero in a finite time, and the disentanglement time is determined by the width of the atomic wave packets. 

This paper is organized as follows. In Sec. II we present an effective Hamiltonian of the total system and then give 
the solution for the eigenequation of the effective Hamiltomian. In Sec. HI we discuss the time evolution of the density 
operator of the total system and obtain the reduced density matrix in the relative-coordinate picture. Furthermore, 
a decoherence factor for arbitrary cavity-field and general atomic spatial states is introduced. In Sec. IV, under the 
assumption that the two atoms' relative position is initially in a superposition state of two Gaussian wave packets, 
we demonstrate the spatial decoherence when the cavity field is initially in a coherent state. In Sec. V we study the 
disentanglement dynamics of the two atoms' internal states. Finally, a concluding summary is given in Sec. VI. 



II. MODEL 



Our system consists of two identical two- level atoms that interact with a single-mode cavity field (Fig. [T]). The two 
atoms are denoted here as atom 1 and atom 2. The mass of each atom is denoted by mo and their atomic transition 
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frequency is wq. Assume tfiat tlie two atoms are spatially separated and located at the positions ii and X2 respectively, 
and the corresponding momenta for atom 1 and atom 2 are pi and p2. Then the momentum for the center of mass 
(cm.) of the two atoms is P = pi +j52, and the relative momentum of the two atoms is p = (pi — p2)/2. The cm. 
position and the relative position are X = {x\ + x^) /2 and x = X\ — X2, respectively. 
Under the rotating-wave approximation, the Hamiltonian of this system reads as 
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) + H.c. 



(1) 



where ct^'^ = \ei) {ei\ - \g^) , cr 
with respect to the excited state 
operators for the cavity field, u and k are the frequency and wave number of the cavity field, and g is the atom-field 
coupling constant. There may be a counter-propagating running- wave mode with a wave vector — fc in such a ring 
cavity, which may affect the atom- field dynamics and contribute to the interaction part. However, as discussed in Ref. 
[Tsl [24 | , we can consider that an atom traverses only an arm of the optical ring cavity and only one cavity mode is 
excited by an external laser. Thus, it is reasonable to consider a single propagating (e.g., a clockwise- running) wave 
mode in a ring cavity. 

We first factorize the evolution operator lJ{t) = cxp{—iHt/h) into a product as 



Ci) {gi\ and irl^' = \gi) (cij (i = 1, 2) are the atomic operators for the i-th atom 
e) and the ground state \g) of the atoms, and a are the creation and annihilation 



U{t) = W{xi)W{x2)Ue{t)W{x2VW{xi)^ 



where W{xi) {i — 1, 2) is a unitary transformation defined by 



W{x, 



exp 



|ei)(ei 



exp 



(2) 



(3) 



which concerns the coupling of the internal levels with the spatial degrees of the atom i. The operator Ue(t) = 
exp (^—iHf,t/H^ is easily proved to be determined by the effective Hamiltonian = Hq + Hi of the system with 
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Here, VLi = ujq + pik/niQ and 1^2 — (^o + P2k/mQ. In the following discussion, we consider that the atomic transition 
is resonant with the cavity mode, i.e., ujq — uj. In this case, we have ujq = ck (c is the light velocity in vacuum). 
Note that pik/niQ = Vik {i = 1,2), where Vi is the velocity of the atom i. Thus, when the velocity of the atom i 




FIG. 1: (Color online) A schematic diagram of the system considered in this paper. A single-mode ring cavity contains two 
atoms 1 and 2. Here, p is the relative momentum of the two atoms and x is the relative position of the two atoms. 
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is far smaller than the light velocity in vacuum (i.e., Vi << c), we have pik/mo << ujq. On the other hand, the 
condition pik/mQ << g needs to be met, such that the influence of the momentum-dependent energy shifts of the 
atomic internal levels on the system dynamics is negligibly small. To see the availability of this condition, let us 
consider a ^^Rb atom with two circular Rydberg levels \g) and |e) (corresponding to principal quantum numbers 50 
and 51). The transition frequency between \g) and |e) is ^ 51.1 GHz |25j. Thus, we have k ^ 10"^ m^^ for the 
case when the transition between \g) and |e) is resonant with the cavity mode. The coupling constant g is on the 
order of 10^ Hz [25l - l27j . For the laser-cooled and optically trapped ^'^Rb atom, we can assume that the atom has a 
measured temperature of T '-^ 180 fiK \2^. A simple calculation shows that the atomic velocity is ~ 23 c.m./s. 
Hence, we have Vk ^ 230 Hz << 10^ Hz, which demonstrates that the approximation of the condition \pik/mQ\ << g 
could be satisfied in practice. The analysis given here shows that the term pik/niQ in fi^ can be neglected, leading to 
Hi = ujq + pik/niQ « Wq. We note that the same approximation was made in Ref. (20| . 

The total excitation number a'I'a-f |ei)(ei| -I- |e2)(e2| is conserved during the interaction. By resolving the eigenequa- 
tion of Hi, in the subspace spanned by states with the total excitation number (n -I- 2), the following eigenstates of 
Hi are obtained: 

1^-)^"^ = \^/2n|ei,e2,n)- V2/i„|5i,52,n + 2), (6) 
1*)^"^ = ^\gi,e2,n + l)-^\ei,g2,n+l), (7) 
1^)3"^ = /in |ei,e2,ri) + i |.gi,e2,n-l- 1) 

+ ^ \ei,g2,n + l) + f2n \gi,g2,n + 2) , (8) 

= -/in |ei, 62, n) i 131,62, n + 1) 

+ ^ |6l,52,?^+ 1) - /2n 131,32,'^ + 2) , (9) 



with the eigenvalues 



e["^ = in + l)tuj = Eo, 4"^=^o + ft^n, e["^ = Eo - hA,,, (10) 



where An = \/%2n -I- 3)(?, /i„ = ^/(n + 1) / [2 (2n -|- 3)], jm = \J in + 2) / [2 (2n -|- 3)] and n is an arbitrary non- 
negative integer. 



III. DECOHERENCE FACTOR 



In this section, we discuss the time evolution of the system, and investigate the spatial decoherence factor for the 
atom-atom relative position. 

Assume that the initial density operator of the whole system is given by 

p(0)=p,(0) ft(0) P/(0), (11) 

where Ps (0) = |'0 (0)) (^/^ (0)| is the initial density operator for the spatial motion of the two atoms, (0) — 
|(/)(0)) ((/'(O)l is the initial density operator for the internal state of the two atoms, and Pj^(O) is the initial den- 
sity operator for the cavity field. 

Let us now assume that the state j?/' (0)) is expanded (in the momentum representation) as 

1^(0))=// dpidp2Cp,p,\pi,p2), (12) 

where \pi) is the momentum eigenstate of atom 1, \p2) is the momentum eigenstate of atom 2, and Cp-^p^ is the 
distribution function satisfying the normalization condition J J^^dpidp2 \Cpj^p^\^ = 1. The state |0(O)) is assumed 
to be a pure separable state |6i,e2). In terms of Fock states \n) and \n') , the density operator (0) is, in general, 
written as 

Pf (0) = E^"-"' I") ("'I ■ (13) 
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The time evolution of p{0) can be obtained by p{t) = U{t)p{0)U^{t). Following the above eigenstates ([B])-®, 
eigenvalues (fTOt of Hi, and based on Eqs. ([2][3]) and (IlllllSp . we find that the density operator p (t) can be written as 

p(t) = ^c„,„,|vl/„(t))(vl/„,(t)|, (14) 

n,n' 



where 

'it 



= / / dpidp2 Cp,p^ exp ■ 



x[Di {n,t) \pi,p2) ® |ei,e2,n) + D2 {n,t) {\pi - hk,p2) ® \gi,e2,n + 1) 

+ \Pi,P2 - hk) (g) \ei,g2,n + 1)) + D3 {n,t) \pi - hk,p2 - hk) ® \gi,g2,n + 2)] (15) 

and 

Di{n,t) = 2/2„cos(A„0 + 2/|„, (16) 
D2{n,t) = -i/i„sin(A„t), (17) 

D3{n,t) = 2/i„/2n[cos(A„t) - 1]. (18) 

We assume that the distribution function Cp-^p^ in Eq. ([T2|) satisfies Cp-^p^ — CpCp] i.e., the initial state |'0 (0)) for 
the spatial motion of the two atoms can be written as IV" (0)) = ® with |^ (0)) = J^^dP Cp\P) 

describing the initial cm. state of the atoms with a momentum distribution function Cp corresponding to the cm. 
momentum eigenstate |P), and |'y3(0)) = J^^dp Cp\p) describing the initial relative position state of the atoms 
with a distribution function Cp corresponding to the relative momentum eigenstate \p). Both Cp and Cp satisfy the 
normalization condition, i.e., |C'p|^ dP — 1 and \Cp\^ dp ^ 1. 

By tracing over the cavity field, the cm. motion, and the internal states of the atoms, we obtain from Eq. (jl4p the 
following elements of the reduced density matrix (in a relative coordinate picture) 

p {x, x', t)^ip (x, t) ip* {x', t) F (x, x', t) , (19) 

where ip {x, t) = dp Cp exp {—ip^t/ {moh) + ipx/h} is the free-evolution state of \(p (0)) expressed in a relative 
coordinate picture, and F{x,x',t) is the decoherence factor, which is given by 

00 

F (x, x\ t) = ^c„,„ \^Dl {n, t) + Dl (n, t)+2 \D2 {n, t)\^ cos [k{x - x')/2]} . (20) 

The decoherence factor F{x,x',t) given here will be used below in our analysis of the spatial decoherence of two 
atomic wave packets. 

IV. SPATIAL DECOHERENCE OF TWO ATOMIC WAVE PACKETS 

We consider that the initial state \(p (0)) of the two atoms' relative position is a superposition of two Gaussian wave 
packets centered at a and —a, respectively; i.e., we have (in the x representation) 

^(a;,0) = -i=[G+(x)+G_(x)], (21) 

— 1/2 

where G±{x) = {V2TTd) ' exp [-{x±af/ (4^^)] and 5 = 1-Hexp [-aV (2^^)]. Here 5 is a normalization constant 
and, for simplicity, we have assumed that the two Gaussian distributions have the same spread d, which is limited in 
dndn < d < a. Here, d^un is much smaller than the wavelength A of the cavity field, but a zero spread is not permitted 
to avoid the atoms indistinguishable by quantum dispersion [29]. 

Assume that the time for the atoms staying in the cavity is so short that we can assume the spatial distribution 
remains unchanged with respect to the initial spatial state for the atomic free evolution. Then Eq. ()19|) can be written 
as a product of the initial state and the decoherence factor 



pix,x',t) ~ ip{x,0) ip*{x',0) F{x,x',t) = p{x,x',0) F{x,x',t) , 



(22) 



5 



where p{x,x',0) = (p {x , 0) cp* {x' , 0) describes the initial state of the atomic relative position. The initial state is 
illustrated in Fig. (3] Here, the two peaks, along the x' = x direction, correspond to the diagonal terms of p{x,x',0), 
while the other two peaks along the x' — —x direction correspond to the off-diagonal terms of p{x,x',0) , which 
represent the coherence between the two wave packets. 

Following Ea.pO|. it can be proved that for the case of x' — x, F {x, x, t) — 1, thus we have p {x, x,t) — p (x, x, 0), 
i.e., the diagonal terms of the density matrix p (x, x' , t) remain unchanged during the time evolution. However, for the 
off-diagonal terms with x' = —x ^ mX {m is an arbitrary nonzero integer), it can be seen that \F {x' = —x,t)\ < 1. 
Thus we will mainly analyze the evolution of the off-diagonal terms of the reduced density matrix below. 

For the cavity field initially in a coherent state |a), we have c„^„ = exp(— |a| 
number. Thus, Eq. (pUj) becomes 



a /n\, where a is a complex 



F(x,x',t) = e-l"l'^^ I Dl {n,t)+2\D2 (n,i)pcos 



k{x — x') 



DIM 



(23) 



Along the x' — —x direction, the decoherence factor takes the following form 



Fix' 



-x,t) 



n=0 



^^Df (n, t)+2 \D2 {n, cos (kx) + Dl (n, t)} , 



(24) 



which is shown in Fig. [31 It can be seen that: 

(i) for X — mA, F{x' — —x,t) remains equal to 1 during the time evolution [Fig. [3ja)]; 

(ii) for X ^ mA, F(x' = —x,t) decays to a constant after a finite time and the value of the constant is determined 
by the value of x [Fig. WL^)] ■ 

(iii) For a given time < , F{x' = —x,t) oscillates with a; in a cosine law [Fig. |3Jc)] , and for the case of a; = {m + 1/2) A, 
the decay of F{x' = —x,t) reaches the maximum with F{x' = —x,gt ^ 1.5) « 0.5 [Fig. [Sjc)]. 

Above, we presented an analysis on the decoherence factor, for the case of the initial state (p{x, 0) described by Eq. 
(PT|). According to the above analysis, it can be concluded that: 

(i) for the case a = mA, the two peaks along the x' — —x direction remain unchanged during the time evolution. 
So the relative motion of the two atoms decouples with the cavity field, and thus there is no decoherence induced by 
the photon recoil. 

(ii) for the case a ^ mX, the two peaks along the x' — ~x direction partially decay in a finite time and there is no 
revival, and the decay of the two peaks depends on a in a cosine law, with a maximum decay at a = (m -I- 1/2) A. 

The evolution of p{x — —x, t) for the initial state •p{x^ 0) with a ^ mA is shown in Fig. Ul which demonstrates 
that the initial state 'p{x^ 0) undergoes a partial decoherence after a finite time and a revival does not occur. Thus, 
the coherence of the atomic relative position is strongly destroyed. The wigner functions for the initial density matrix 
p{x, X , 0) shown in Fig. [2] and the density matrix p{x, x , t) given in Eq. (|22|) are shown in Fig. [5j Strong oscillations 
together with negative values [Fig. [S^a)] indicate quantum coherence between the two wave packets, and it turns 




FIG. 2: (Color online) Plot of the density matrix p{x,x',0) representing the initial state for the relative position of the two 
atoms. It is a superposition of two Gaussian wave packets. The two peaks along the x' — ~x diagonal direction represent the 
coherence between the two wave packets. 
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FIG. 3: (Color online) (a) Evolution of the decoherence factor F{x' = ~x, t) with time t and relative position x. Decoherence 
factor F{x' — —x,t) for (b) some selected relative positions and (c) a selected time, (a), (b) and (c) were ploted for the cavity 
field initially in a coherent state with a = 10. 

out that the oscillations are partially damped by decoherence [Fig. Efb)]. This may be induced by the entanglement 
between the two-atom spatial motion and the cavity field resulting from the photonic back-action. This entangling 
process inevitably destroys the coherence of the spatial motion of the two atoms. Since the single-mode cavity field 
is just an environment with a few degrees of freedom, the spatial coherence is not completely destroyed. 



V. DISENTANGLEMENT DYNAMICS OF TWO ATOMS 



To begin with, let us assume that: 

(i) the initial spatial state \ip (0)) of the two atoms is a separable state 

|V(0)) = |m)i®|m)2- (25) 
Here, \fi)-^ and |/i)2 represent two Gaussian wave packets describing the spatial distribution of the two atoms, respec- 



FIG. 4: (Color online) Density matrix p{x' — —x,t), which evolves from the initial state (|2ip with a 7^ mX for the cavity field 
initially in a coherent state with a — 10. 




FIG. 5: (Color online) (a) The Wigner function for the initial density matrix p{x,x',0) shown in Fig. [2] (b) The Wigner 
function for the density matrix in Eq. (|22p whose off-diagonal terms were shown in Fig. [l] (a) and (b) were plotted for gt = 2. 
Here, p is the relative momentum of the two atoms and x is the relative position of the two atoms. 



tively. In the coordinate representation, \fi)^ is expressed as (i = 1,2): 



dpi Cp^ exp ( -^PiXi 



2t:(P 



1/4 



exp 



{xi + aif 



(26) 



where is the center of the Gaussian function /i (2:^,0) and d is the width of the Gaussian function /i {xi^O). The 
coefficient Cp^ is given by 



— 



exp 



(27) 



For convenience, we assume that oi = —02 = a/2. The initial density operator for the spatial motion of the two 
atoms is p, (0) = \^ (0)) (0)|. 
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FIG. 6: (Color online) Time evolution of the concurrence C for the case when the initial internal state of the two atoms is an 
entangled state (0)) = cos 7 Ifli, 52) + sin 7 |ei, 62), with 7 = 7r/4. The dotted line is for e"" with d = a/lQO. 



(ii) the initial internal state (0)) of the two atoms is an entangled state, which is given by 

10(0)) = cos7|gi,g2) +sin7|ei,e2) . 



(28) 



The initial density operator for the internal state of the two atoms is p^ (0) = |0 (0)) (0 (0)| . 
(ill) the cavity field is initially in a vacuum state, i.e., pj (0) — |0) (0|. 

The initial density operator of the whole system is /o(0) = Ps(0) Pi(0) PfW- time t, the density operator 
of the system is p{t) = U{t)p{0) U^{t). Here, U{t) is the unitary operator in Eq. ([2]). After tracing p{t) over the 
electromagnetic field and the spatial degrees of freedom of the atoms, we obtain the following reduced density matrix 
for the internal state of the two atoms: 



P^it) = 



I a w 

6 z 

z* c 

\w* d 



(29) 



which is written in a basis formed by |ei,e2), |ei,(72), 151562) and |5i,52)- Here, 

h = |i?2(0,0)|'sin2 7, 
d = cos2 7 + i:>^(0,0)sin^7, 
w = i:'i(0,0)cos7sin7e-*2'^*e-"(*\ 
z = 6e-«(*), 

with ^(t) = s{€) + d^fc^ — ika and sif) = h^k'^t^ / {Ad'^rn^. The parameter s{t) is obtained after tracing p (t) over the 
atomic spatial degrees of freedom, which causes the strong time dependence of the concurrence as shown below. 
A popular measure of entanglement for p^{t) is given by the concurrence [s^], which is defined as 




where the Xi are the eigenvalues of the non-Hermitian matrix p^^'p^ in descending order, and p^ = 
{p'y®'^y)Pi {.^v®^v) ■ Here, p* is the complex conjugate of p.^, and tr^ is the usual Pauli operator for atom i 
{i — 1,2). The concurrence C varies from C = for a non-entangled state to C = 1 for a maximally-entangled state. 
The concurrence of the density matrix p^ (t) in Eq. is [sil . [33 ] 

C {p,{t)) = 2 max{0, \z\ - Vad, \w\ - Vfoc}, (31) 

which is shown in Fig. |6l Eq. tells us that the concurrence evolves with an non-smoothly exponential decay due 
to the factor as shown in Fig. [H 
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FIG. 7: (Color online) Time evolution of the concurrence C for the case when the initial internal state of the two atoms is an 
entangled state (0)) — cos 7 lei, ^2) + sin 7 jgri, 62), with 7 = 7r/4. The dotted line is for e"" with d = a/100. 



Alternatively, the initial internal state for the two atoms could be another type of entangled state given by 

10(0)) = cos7|ei,52) +sin7|gi,e2) . (32) 

By doing a calculation similar to the one above, it can be seen that for the entanglement state |0 (0)) considered here, 
the reduced density matrix p^{t) for the internal state of the two atoms is 

/ \ 
\0 wj 

where 

b = Al{t) cos^ -/ + A^_{t) si-n^ J 

+2 {t)A- (t) cos 7 sin je'^ cos (ak) , 
c = Aiit) cos'^ J + A'^_^_(t) sin^ J 

+2 {t)A^ (t) cos 7 sin 76""^ cos(afc) , 
w = |i3(t)p + 2cos7sin7e~''cos(afc) , 
z = e~''^''>{A+{t)A^{t)[cos^je-°'+sm^je-^] 

+ [Alit) + A^(t)e-''] cos 7 sin 7} , 

with 

cos(V2.9^) ± 1 
A±{t) ^ , 

B{t) ^ -ism{V2gt)/V2. 

Here, S = (Pk'^, a = (Pk"^ + iak, (i = d^fc^ — iak and 77 = AtPk"^ + i2ak. 

The concurrence for the reduced density matrix p^it) in Eq. is given by C {p^{t)) = 2 |z| , which is shown in 
Fig. [71 One can see that there is a common factor e"**^*-' with s{t) = h?k'^t^/ (4ci^mg), appearing in the expression of 
z above. Because of this factor e~^^^\ the value of \z\ decays to zero with time asymptotically. Here, Fig. Eland Fig. 
[71 show that the oscillations of the entanglement are damped exponentially, but the decay of the entanglement occurs 
faster when the width d of the initial Gaussian function becomes smaller, and the entanglement decays to almost 
zero in a finite time. So the back-action of the emitted photons may induce an entanglement between the two-atom 
internal states and the two-atom spatial motion states, and this correlation destroys the coherence of the internal 
states of the two atoms. The smaller the width of the wave packets, the larger the uncertainty of the atomic spatial 
momentum, which results in a faster destruction of the atomic internal-state entanglement. 
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VI. CONCLUSION 

We have considered a system of two atoms interacting with a single-mode cavity field. Since the spatial degrees 
of freedom of the atoms arc considered, the Hamiltonian of the whole system becomes comphcated. To solve the 
Schordinger equation, we have introduced two unitary transformations (involving the coupling of the internal levels 
with the spacial degrees of the atoms), presented an effective Hamiltonian of the system, and given the analytic 
solutions to the eigenequation of the effective Hamiltonian. Based on these, we have investigated the decoherence 
dynamics for the relative position of the two atoms, and presented a decoherence factor for a general cavity-field state 
and an arbitrary atomic spatial state, which might be useful for future related works. 

Under the assumption that the atomic relative position is in a superposition of two Gaussian wave packets, we 
have demonstrated the spatial decoherence of the two atoms' relative position, for a cavity field initially in a coherent 
state. Our results show that the spatial decoherence of the two atoms depends strongly on the relative position of 
the two atoms. Interestingly, we found that when the relative position of the two atoms is an integral multiple of the 
wavelength of the cavity field, the spatial coherence of the relative position of the two atoms is not destroyed by the 
photon recoil; However, when the relative position of the two atoms is not an integral multiple of the wavelength of 
the cavity field, spatial decoherence of the relative position of the two atoms happens. 

Furthermore, we have studied the entanglement dynamics of the internal states of the two atoms interacting with 
a single-mode cavity field. Our results show that the entanglement, measured by the concurrence, decays to almost 
zero in a finite time. Thus, the back-action of the photons emitted from the two atoms may be a fundamental process 
destroying the cintanglement of atoms. 

From this work, it can be concluded that there exists a phenomenon that the moving qubits (e.g., atoms) placed in 
a cavity, may suffer from the spatial decoherence. Therefore, it is important to overcome the influence of the spatial 
motion of qubits on their entanglement. 
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